Analysis of the potential relationship between COVID-19 and Behcet's disease using transcriptome data.

To investigate the potential role of COVID-19 in relation to Behcet's disease (BD) and to search for relevant biomarkers. We used a bioinformatics approach to download transcriptomic data from peripheral blood mononuclear cells (PBMCs) of COVID-19 patients and PBMCs of BD patients, screened the common differential genes between COVID-19 and BD, performed gene ontology (GO) and pathway analysis, and constructed the protein-protein interaction (PPI) network, screened the hub genes and performed co-expression analysis. In addition, we constructed the genes-transcription factors (TFs)-miRNAs network, the genes-diseases network and the genes-drugs network to gain insight into the interactions between the 2 diseases. We used the RNA-seq dataset from the GEO database (GSE152418, GSE198533). We used cross-analysis to obtain 461 up-regulated common differential genes and 509 down-regulated common differential genes, mapped the PPI network, and used Cytohubba to identify the 15 most strongly associated genes as hub genes (ACTB, BRCA1, RHOA, CCNB1, ASPM, CCNA2, TOP2A, PCNA, AURKA, KIF20A, MAD2L1, MCM4, BUB1, RFC4, and CENPE). We screened for statistically significant hub genes and found that ACTB was in low expression of both BD and COVID-19, and ASPM, CCNA2, CCNB1, and CENPE were in low expression of BD and high expression of COVID-19. GO analysis and pathway analysis was then performed to obtain common pathways and biological response processes, which suggested a common association between BD and COVID-19. The genes-TFs-miRNAs network, genes-diseases network and genes-drugs network also play important roles in the interaction between the 2 diseases. Interaction between COVID-19 and BD exists. ACTB, ASPM, CCNA2, CCNB1, and CENPE as potential biomarkers for 2 diseases.

Exploration of the potential common pathogenic mechanisms in COVID-19 and silicosis by using bioinformatics and system biology.

Silicosis is an occupational lung disease that is common worldwide. In recent years, coronavirus disease 2019 (COVID-19) has provided daunting challenges to public healthcare systems globally. Although multiple studies have shown a close link between COVID-19 and other respiratory diseases, the inter-relational mechanisms between COVID-19 and silicosis remain unclear. This study aimed to explore the shared molecular mechanisms and drug targets of COVID-19 and silicosis. Gene expression profiling identified four modules that were most closely associated with both diseases. Furthermore, we performed functional analysis and constructed a protein-protein interaction network. Seven hub genes (budding uninhibited by benzimidazoles 1 [BUB1], protein regulator of cytokinesis 1 [PRC1], kinesin family member C1 [KIFC1], ribonucleotide reductase regulatory subunit M2 [RRM2], cyclin-dependent kinase inhibitor 3 [CDKN3], Cyclin B2 [CCNB2], and minichromosome maintenance complex component 6 [MCM6]) were involved in the interaction between COVID-19 and silicosis. We investigated how diverse microRNAs and transcription factors regulate these seven genes. Subsequently, the correlation between the hub genes and infiltrating immune cells was explored. Further in-depth analyses were performed based on single-cell transcriptomic data from COVID-19, and the expression of hub-shared genes was characterized and located in multiple cell clusters. Finally, molecular docking results reveal small molecular compounds that may improve COVID-19 and silicosis. The current study reveals the common pathogenesis of COVID-19 and silicosis, which may provide a novel reference for further research.

Cell-attribute aware community detection improves differential abundance testing from single-cell RNA-Seq data.

Variations of cell-type proportions within tissues could be informative of biological aging and disease risk. Single-cell RNA-sequencing offers the opportunity to detect such differential abundance patterns, yet this task can be statistically challenging due to the noise in single-cell data, inter-sample variability and because such patterns are often of small effect size. Here we present a differential abundance testing paradigm called ELVAR that uses cell attribute aware clustering when inferring differentially enriched communities within the single-cell manifold. Using simulated and real single-cell and single-nucleus RNA-Seq datasets, we benchmark ELVAR against an analogous algorithm that uses Louvain for clustering, as well as local neighborhood-based methods, demonstrating that ELVAR improves the sensitivity to detect cell-type composition shifts in relation to aging, precancerous states and Covid-19 phenotypes. In effect, leveraging cell attribute information when inferring cell communities can denoise single-cell data, avoid the need for batch correction and help retrieve more robust cell states for subsequent differential abundance testing. ELVAR is available as an open-source R-package.

Bioinformatics and system biology approach to identify the influences among COVID-19, ARDS and sepsis.

Background Severe coronavirus disease 2019 (COVID -19) has led to severe pneumonia or acute respiratory distress syndrome (ARDS) worldwide. we have noted that many critically ill patients with COVID-19 present with typical sepsis-related clinical manifestations, including multiple organ dysfunction syndrome, coagulopathy, and septic shock. The molecular mechanisms that underlie COVID-19, ARDS and sepsis are not well understood. The objectives of this study were to analyze potential molecular mechanisms and identify potential drugs for the treatment of COVID-19, ARDS and sepsis using bioinformatics and a systems biology approach. Methods Three RNA-seq datasets (GSE171110, GSE76293 and GSE137342) from Gene Expression Omnibus (GEO) were employed to detect mutual differentially expressed genes (DEGs) for the patients with the COVID-19, ARDS and sepsis for functional enrichment, pathway analysis, and candidate drugs analysis. Results We obtained 110 common DEGs among COVID-19, ARDS and sepsis. ARG1, FCGR1A, MPO, and TLR5 are the most influential hub genes. The infection and immune-related pathways and functions are the main pathways and molecular functions of these three diseases. FOXC1, YY1, GATA2, FOXL, STAT1 and STAT3 are important TFs for COVID-19. mir-335-5p, miR-335-5p and hsa-mir-26a-5p were associated with COVID-19. Finally, the hub genes retrieved from the DSigDB database indicate multiple drug molecules and drug-targets interaction. Conclusion We performed a functional analysis under ontology terms and pathway analysis and found some common associations among COVID-19, ARDS and sepsis. Transcription factors-genes interaction, protein-drug interactions, and DEGs-miRNAs coregulatory network with common DEGs were also identified on the datasets. We believe that the candidate drugs obtained in this study may contribute to the effective treatment of COVID-19.

Using bioinformatics and systems biology methods to identify the mechanism of interaction between COVID-19 and nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is considered a risk factor for severe COVID-19, but the mechanism remains unknown. This study used bioinformatics to help define the relationship between these diseases. The GSE147507 (COVID-19), GSE126848 (NAFLD), and GSE63067 (NAFLD-2) datasets were screened using the Gene Expression Omnibus. Common differentially expressed genes were then identified using a Venn diagram. Gene ontology analysis and KEGG pathway enrichment were performed on the differentially expressed genes. A protein-protein interaction network was also constructed using the STRING platform, and key genes were identified using the Cytoscape plugin. GES63067 was selected for validation of the results. Analysis of ferroptosis gene expression during the development of the 2 diseases and prediction of their upstream miRNAs and lncRNAs. In addition, transcription factors (TFs) and miRNAs related to key genes were identified. Effective drugs that act on target genes were found in the DSigDB. The GSE147507 and GSE126848 datasets were crossed to obtain 28 co-regulated genes, 22 gene ontology terms, 3 KEGG pathways, and 10 key genes. NAFLD may affect COVID-19 progression through immune function and inflammatory signaling pathways. CYBB was predicted to be a differential ferroptosis gene associated with 2 diseases, and the CYBB-hsa-miR-196a/b-5p-TUG1 regulatory axis was identified. TF-gene interactions and TF-miRNA coregulatory network were constructed successfully. A total of 10 drugs, (such as Eckol, sulfinpyrazone, and phenylbutazone) were considered as target drugs for Patients with COVID-19 and NAFLD. This study identified key gene and defined molecular mechanisms associated with the progression of COVID-19 and NAFLD. COVID-19 and NAFLD progression may regulate ferroptosis through the CYBB-hsa-miR-196a/b-5p-TUG1 axis. This study provides additional drug options for the treatment of COVID-19 combined with NAFLD disease.

Bioinformatics and system biology approach to identify potential common pathogenesis for COVID-19 infection and osteoarthritis.

A growing of evidence has showed that patients with osteoarthritis (OA) had a higher coronavirus 2019 (COVID-19) infection rate and a poorer prognosis after infected it. Additionally, scientists have also discovered that COVID-19 infection might cause pathological changes in the musculoskeletal system. However, its mechanism is still not fully elucidated. This study aims to further explore the sharing pathogenesis of patients with both OA and COVID-19 infection and find candidate drugs. Gene expression profiles of OA (GSE51588) and COVID-19 (GSE147507) were obtained from the Gene Expression Omnibus (GEO) database. The common differentially expressed genes (DEGs) for both OA and COVID-19 were identified and several hub genes were extracted from them. Then gene and pathway enrichment analysis of the DEGs were performed; protein-protein interaction (PPI) network, transcription factor (TF)-gene regulatory network, TF-miRNA regulatory network and gene-disease association network were constructed based on the DEGs and hub genes. Finally, we predicted several candidate molecular drugs related to hub genes using DSigDB database. The receiver operating characteristic curve (ROC) was applied to evaluate the accuracy of hub genes in the diagnosis of both OA and COVID-19. In total, 83 overlapping DEGs were identified and selected for subsequent analyses. CXCR4, EGR2, ENO1, FASN, GATA6, HIST1H3H, HIST1H4H, HIST1H4I, HIST1H4K, MTHFD2, PDK1, TUBA4A, TUBB1 and TUBB3 were screened out as hub genes, and some showed preferable values as diagnostic markers for both OA and COVID-19. Several candidate molecular drugs, which are related to the hug genes, were identified. These sharing pathways and hub genes may provide new ideas for further mechanistic studies and guide more individual-based effective treatments for OA patients with COVID-19 infection.

Facilitating systems-level analyses of all-cause and Covid-mediated sepsis through SeptiSearch, a manually-curated compendium of dysregulated gene sets.

Background: Sepsis is a dysfunctional host response to infection. The syndrome leads to millions of deaths annually (19.7% of all deaths in 2017) and is the cause of most deaths from severe Covid infections. High throughput sequencing or 'omics' experiments in molecular and clinical sepsis research have been widely utilized to identify new diagnostics and therapies. Transcriptomics, quantifying gene expression, has dominated these studies, due to the efficiency of measuring gene expression in tissues and the technical accuracy of technologies like RNA-Seq. Objective: Most of these studies seek to uncover novel mechanistic insights into sepsis pathogenesis and diagnostic gene signatures by identifying genes differentially expressed between two or more relevant conditions. However, little effort has been made, to date, to aggregate this knowledge from such studies. In this study we sought to build a compendium of previously described gene sets that combines knowledge gained from sepsis-associated studies. This would enable the identification of genes most associated with sepsis pathogenesis, and the description of the molecular pathways commonly associated with sepsis. Methods: PubMed was searched for studies using transcriptomics to characterize acute infection/sepsis and severe sepsis (i.e., sepsis combined with organ failure). Several studies were identified that used transcriptomics to identify differentially expressed (DE) genes, predictive/prognostic signatures, and underlying molecular responses and pathways. The molecules included in each gene set were collected, in addition to the relevant study metadata (e.g., patient groups used for comparison, sample collection time point, tissue type, etc.). Results: After performing extensive literature curation of 74 sepsis-related publications involving transcriptomics, 103 unique gene sets (comprising 20,899 unique genes) from thousands of patients were collated together with associated metadata. Frequently described genes included in gene sets as well as the molecular mechanisms they were involved in were identified. These mechanisms included neutrophil degranulation, generation of second messenger molecules, IL-4 and -13 signaling, and IL-10 signaling among many others. The database, which we named SeptiSearch, is made available in a web application created using the Shiny framework in R, (available at https://septisearch.ca). Conclusions: SeptiSearch provides members of the sepsis community the bioinformatic tools needed to leverage and explore the gene sets contained in the database. This will allow the gene sets to be further scrutinized and analyzed for their enrichment in user-submitted gene expression data and used for validation of in-house gene sets/signatures.

Integrative Analysis of Metabolomic and Transcriptomic Data Reveals the Mechanism of Color Formation in Corms of Pinellia ternata.

Pinellia ternata (Thunb.) Breit. (P. ternata) is a very important plant that is commonly used in traditional Chinese medicine. Its corms can be used as medicine and function to alleviate cough, headache, and phlegm. The epidermis of P. ternata corms is often light yellow to yellow in color; however, within the range of P. ternata found in JingZhou City in Hubei Province, China, there is a form of P. ternata in which the epidermis of the corm is red. We found that the total flavonoid content of red P. ternata corms is significantly higher than that of yellow P. ternata corms. The objective of this study was to understand the molecular mechanisms behind the difference in epidermal color between the two forms of P. ternata. The results showed that a high content of anthocyanidin was responsible for the red epidermal color in P. ternata, and 15 metabolites, including cyanidin-3-O-rutinoside-5-O-glucoside, cyanidin-3-O-glucoside, and cyanidin-3-O-rutinoside, were screened as potential color markers in P. ternata through metabolomic analysis. Based on an analysis of the transcriptome, seven genes, including PtCHS1, PtCHS2, PtCHI1, PtDFR5, PtANS, PtUPD-GT2, and PtUPD-GT3, were found to have important effects on the biosynthesis of anthocyanins in the P. ternata corm epidermis. Furthermore, two transcription factors (TFs), bHLH1 and bHLH2, may have regulatory functions in the biosynthesis of anthocyanins in red P. ternata corms. Using an integrative analysis of the metabolomic and transcriptomic data, we identified five genes, PtCHI, PtDFR2, PtUPD-GT1, PtUPD-GT2, and PtUPD-GT3, that may play important roles in the presence of the red epidermis color in P. ternata corms.

SARS-CoV-2 airway infection results in the development of somatosensory abnormalities in a hamster model.

Although largely confined to the airways, SARS-CoV-2 infection has been associated with sensory abnormalities that manifest in both acute and chronic phenotypes. To gain insight on the molecular basis of these sensory abnormalities, we used the golden hamster model to characterize and compare the effects of infection with SARS-CoV-2 and influenza A virus (IAV) on the sensory nervous system. We detected SARS-CoV-2 transcripts but no infectious material in the cervical and thoracic spinal cord and dorsal root ganglia (DRGs) within the first 24 hours of intranasal virus infection. SARS-CoV-2-infected hamsters exhibited mechanical hypersensitivity that was milder but prolonged compared with that observed in IAV-infected hamsters. RNA sequencing analysis of thoracic DRGs 1 to 4 days after infection suggested perturbations in predominantly neuronal signaling in SARS-CoV-2-infected animals as opposed to type I interferon signaling in IAV-infected animals. Later, 31 days after infection, a neuropathic transcriptome emerged in thoracic DRGs from SARS-CoV-2-infected animals, which coincided with SARS-CoV-2-specific mechanical hypersensitivity. These data revealed potential targets for pain management, including the RNA binding protein ILF3, which was validated in murine pain models. This work elucidates transcriptomic signatures in the DRGs triggered by SARS-CoV-2 that may underlie both short- and long-term sensory abnormalities.

Identifying key genes related to inflammasome in severe COVID-19 patients based on a joint model with random forest and artificial neural network.

Background: The coronavirus disease 2019 (COVID-19) has been spreading astonishingly and caused catastrophic losses worldwide. The high mortality of severe COVID-19 patients is an serious problem that needs to be solved urgently. However, the biomarkers and fundamental pathological mechanisms of severe COVID-19 are poorly understood. The aims of this study was to explore key genes related to inflammasome in severe COVID-19 and their potential molecular mechanisms using random forest and artificial neural network modeling. Methods: Differentially expressed genes (DEGs) in severe COVID-19 were screened from GSE151764 and GSE183533 via comprehensive transcriptome Meta-analysis. Protein-protein interaction (PPI) networks and functional analyses were conducted to identify molecular mechanisms related to DEGs or DEGs associated with inflammasome (IADEGs), respectively. Five the most important IADEGs in severe COVID-19 were explored using random forest. Then, we put these five IADEGs into an artificial neural network to construct a novel diagnostic model for severe COVID-19 and verified its diagnostic efficacy in GSE205099. Results: Using combining P value < 0.05, we obtained 192 DEGs, 40 of which are IADEGs. The GO enrichment analysis results indicated that 192 DEGs were mainly involved in T cell activation, MHC protein complex and immune receptor activity. The KEGG enrichment analysis results indicated that 192 GEGs were mainly involved in Th17 cell differentiation, IL-17 signaling pathway, mTOR signaling pathway and NOD-like receptor signaling pathway. In addition, the top GO terms of 40 IADEGs were involved in T cell activation, immune response-activating signal transduction, external side of plasma membrane and phosphatase binding. The KEGG enrichment analysis results indicated that IADEGs were mainly involved in FoxO signaling pathway, Toll-like receptor, JAK-STAT signaling pathway and Apoptosis. Then, five important IADEGs (AXL, MKI67, CDKN3, BCL2 and PTGS2) for severe COVID-19 were screened by random forest analysis. By building an artificial neural network model, we found that the AUC values of 5 important IADEGs were 0.972 and 0.844 in the train group (GSE151764 and GSE183533) and test group (GSE205099), respectively. Conclusion: The five genes related to inflammasome, including AXL, MKI67, CDKN3, BCL2 and PTGS2, are important for severe COVID-19 patients, and these molecules are related to the activation of NLRP3 inflammasome. Furthermore, AXL, MKI67, CDKN3, BCL2 and PTGS2 as a marker combination could be used as potential markers to identify severe COVID-19 patients.

Exploring Potential Biomarkers and Molecular Mechanisms of Ischemic Cardiomyopathy and COVID-19 Comorbidity Based on Bioinformatics and Systems Biology.

Cardiovascular complications combined with COVID-19 (SARS-CoV-2) lead to a poor prognosis in patients. The common pathogenesis of ischemic cardiomyopathy (ICM) and COVID-19 is still unclear. Here, we explored potential molecular mechanisms and biomarkers for ICM and COVID-19. Common differentially expressed genes (DEGs) of ICM (GSE5406) and COVID-19 (GSE164805) were identified using GEO2R. We performed enrichment and protein-protein interaction analyses and screened key genes. To confirm the diagnostic performance for these hub genes, we used external datasets (GSE116250 and GSE211979) and plotted ROC curves. Transcription factor and microRNA regulatory networks were constructed for the validated hub genes. Finally, drug prediction and molecular docking validation were performed using cMAP. We identified 81 common DEGs, many of which were enriched in terms of their relation to angiogenesis. Three DEGs were identified as key hub genes (HSP90AA1, HSPA9, and SRSF1) in the protein-protein interaction analysis. These hub genes had high diagnostic performance in the four datasets (AUC > 0.7). Mir-16-5p and KLF9 transcription factor co-regulated these hub genes. The drugs vindesine and ON-01910 showed good binding performance to the hub genes. We identified HSP90AA1, HSPA9, and SRSF1 as markers for the co-pathogenesis of ICM and COVID-19, and showed that co-pathogenesis of ICM and COVID-19 may be related to angiogenesis. Vindesine and ON-01910 were predicted as potential therapeutic agents. Our findings will contribute to a deeper understanding of the comorbidity of ICM with COVID-19.

Gene expression in organoids: an expanding horizon.

Recent development of human three-dimensional organoid cultures has opened new doors and opportunities ranging from modelling human development in vitro to personalised cancer therapies. These new in vitro systems are opening new horizons to the classic understanding of human development and disease. However, the complexity and heterogeneity of these models requires cutting-edge techniques to capture and trace global changes in gene expression to enable identification of key players and uncover the underlying molecular mechanisms. Rapid development of sequencing approaches made possible global transcriptome analyses and epigenetic profiling. Despite challenges in organoid culture and handling, these techniques are now being adapted to embrace organoids derived from a wide range of human tissues. Here, we review current state-of-the-art multi-omics technologies, such as single-cell transcriptomics and chromatin accessibility assays, employed to study organoids as a model for development and a platform for precision medicine.

Reconstruction of the cell pseudo-space from single-cell RNA sequencing data with scSpace.

Tissues are highly complicated with spatial heterogeneity in gene expression. However, the cutting-edge single-cell RNA-seq technology eliminates the spatial information of individual cells, which contributes to the characterization of cell identities. Herein, we propose single-cell spatial position associated co-embeddings (scSpace), an integrative method to identify spatially variable cell subpopulations by reconstructing cells onto a pseudo-space with spatial transcriptome references (Visium, STARmap, Slide-seq, etc.). We benchmark scSpace with both simulated and biological datasets, and demonstrate that scSpace can accurately and robustly identify spatially variated cell subpopulations. When employed to reconstruct the spatial architectures of complex tissue such as the brain cortex, the small intestinal villus, the liver lobule, the kidney, the embryonic heart, and others, scSpace shows promising performance on revealing the pairwise cellular spatial association within single-cell data. The application of scSpace in melanoma and COVID-19 exhibits a broad prospect in the discovery of spatial therapeutic markers.

Identification of biomarkers, pathways, and potential therapeutic targets for heart failure using next-generation sequencing data and bioinformatics analysis.

BACKGROUND: Heart failure (HF) is the most common cardiovascular diseases and the leading cause of cardiovascular diseases related deaths. Increasing molecular targets have been discovered for HF prognosis and therapy. However, there is still an urgent need to identify novel biomarkers. Therefore, we evaluated biomarkers that might aid the diagnosis and treatment of HF. METHODS: We searched next-generation sequencing (NGS) dataset (GSE161472) and identified differentially expressed genes (DEGs) by comparing 47 HF samples and 37 normal control samples using limma in R package. Gene ontology (GO) and pathway enrichment analyses of the DEGs were performed using the g: Profiler database. The protein-protein interaction (PPI) network was plotted with Human Integrated Protein-Protein Interaction rEference (HiPPIE) and visualized using Cytoscape. Module analysis of the PPI network was done using PEWCC1. Then, miRNA-hub gene regulatory network and TF-hub gene regulatory network were constructed by Cytoscape software. Finally, we performed receiver operating characteristic (ROC) curve analysis to predict the diagnostic effectiveness of the hub genes. RESULTS: A total of 930 DEGs, 464 upregulated genes and 466 downregulated genes, were identified in HF. GO and REACTOME pathway enrichment results showed that DEGs mainly enriched in localization, small molecule metabolic process, SARS-CoV infections, and the citric acid tricarboxylic acid (TCA) cycle and respiratory electron transport. After combining the results of the PPI network miRNA-hub gene regulatory network and TF-hub gene regulatory network, 10 hub genes were selected, including heat shock protein 90 alpha family class A member 1 (HSP90AA1), arrestin beta 2 (ARRB2), myosin heavy chain 9 (MYH9), heat shock protein 90 alpha family class B member 1 (HSP90AB1), filamin A (FLNA), epidermal growth factor receptor (EGFR), phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1), cullin 4A (CUL4A), YEATS domain containing 4 (YEATS4), and lysine acetyltransferase 2B (KAT2B). CONCLUSIONS: This discovery-driven study might be useful to provide a novel insight into the diagnosis and treatment of HF. However, more experiments are needed in the future to investigate the functional roles of these genes in HF.

Benchmarking of analytical combinations for COVID-19 outcome prediction using single-cell RNA sequencing data.

The advances of single-cell transcriptomic technologies have led to increasing use of single-cell RNA sequencing (scRNA-seq) data in large-scale patient cohort studies. The resulting high-dimensional data can be summarized and incorporated into patient outcome prediction models in several ways; however, there is a pressing need to understand the impact of analytical decisions on such model quality. In this study, we evaluate the impact of analytical choices on model choices, ensemble learning strategies and integrate approaches on patient outcome prediction using five scRNA-seq COVID-19 datasets. First, we examine the difference in performance between using single-view feature space versus multi-view feature space. Next, we survey multiple learning platforms from classical machine learning to modern deep learning methods. Lastly, we compare different integration approaches when combining datasets is necessary. Through benchmarking such analytical combinations, our study highlights the power of ensemble learning, consistency among different learning methods and robustness to dataset normalization when using multiple datasets as the model input.

Transcriptomic and antiviral analyses of PoIFN-Delta5 against porcine enteric viruses in porcine intestinal epithelial cells.

The interferon-delta family was first reported in domestic pigs and belongs to the type I interferon (IFN-I) family. The enteric viruses could cause diarrhea in newborn piglets with high morbidity and mortality. We researched the function of the porcine IFN-delta (PoIFN-δ) family in the porcine intestinal epithelial cells (IPEC-J2) cells infected with porcine epidemic diarrhea virus (PEDV). Our study found that all PoIFN-δs shared a typical IFN-I signature and could be divided into five branches in the phylogenic tree. Different strains of PEDV could induce typical IFN transitorily, and the virulent strain AH2012/12 had the strongest induction of porcine IFN-δ and IFN-alpha (PoIFN-α) in the early stage of infection. In addition, it was found that PoIFN-δ5/6/9/11 and PoIFN-δ1/2 were highly expressed in the intestine. PoIFN-δ5 had a better antiviral effect on PEDV compared to PoIFN-δ1 due to its higher induction of ISGs. PoIFN-δ1 and PoIFN-δ5 also activated JAK-STAT and IRS signaling. For other enteric viruses, transmissible gastroenteritis virus (TGEV), porcine deltacoronavirus (PDCoV), and porcine rotavirus (PoRV), PoIFN-δ1 and PoIFN-δ5 both showed an excellent antiviral effect. Transcriptome analyses uncovered the differences in host responses to PoIFN-α and PoIFN-δ5 and revealed thousands of differentially expressed genes were mainly enriched in the inflammatory response, antigen processing and presentation, and other immune-related pathways. PoIFN-δ5 would be a potential antiviral drug, especially against porcine enteric viruses. These studies were the first to report the antiviral function against porcine enteric viruses and broaden the new acquaintances of this type of interferon though not novelly discovered.

Bioinformatics Analysis of Candidate Genes and Pathways Related to Hepatocellular Carcinoma in China: A Study Based on Public Databases.

Background and Objective: Hepatocellular carcinoma (HCC) is a highly aggressive malignant tumor of the digestive system worldwide. Chronic hepatitis B virus (HBV) infection and aflatoxin exposure are predominant causes of HCC in China, whereas hepatitis C virus (HCV) infection and alcohol intake are likely the main risk factors in other countries. It is an unmet need to recognize the underlying molecular mechanisms of HCC in China. Methods: In this study, microarray datasets (GSE84005, GSE84402, GSE101685, and GSE115018) derived from Gene Expression Omnibus (GEO) database were analyzed to obtain the common differentially expressed genes (DEGs) by R software. Moreover, the gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed by using Database for Annotation, Visualization and Integrated Discovery (DAVID). Furthermore, the protein-protein interaction (PPI) network was constructed, and hub genes were identified by the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape, respectively. The hub genes were verified using Gene Expression Profiling Interactive Analysis (GEPIA), UALCAN, and Kaplan-Meier Plotter online databases were performed on the TCGA HCC dataset. Moreover, the Human Protein Atlas (HPA) database was used to verify candidate genes' protein expression levels. Results: A total of 293 common DEGs were screened, including 103 up-regulated genes and 190 down-regulated genes. Moreover, GO analysis implied that common DEGs were mainly involved in the oxidation-reduction process, cytosol, and protein binding. KEGG pathway enrichment analysis presented that common DEGs were mainly enriched in metabolic pathways, complement and coagulation cascades, cell cycle, p53 signaling pathway, and tryptophan metabolism. In the PPI network, three subnetworks with high scores were detected using the Molecular Complex Detection (MCODE) plugin. The top 10 hub genes identified were CDK1, CCNB1, AURKA, CCNA2, KIF11, BUB1B, TOP2A, TPX2, HMMR and CDC45. The other public databases confirmed that high expression of the aforementioned genes related to poor overall survival among patients with HCC. Conclusion: This study primarily identified candidate genes and pathways involved in the underlying mechanisms of Chinese HCC, which is supposed to provide new targets for the diagnosis and treatment of HCC in China.

HGCA2.0: An RNA-Seq Based Webtool for Gene Coexpression Analysis in Homo sapiens.

Genes with similar expression patterns in a set of diverse samples may be considered coexpressed. Human Gene Coexpression Analysis 2.0 (HGCA2.0) is a webtool which studies the global coexpression landscape of human genes. The website is based on the hierarchical clustering of 55,431 Homo sapiens genes based on a large-scale coexpression analysis of 3500 GTEx bulk RNA-Seq samples of healthy individuals, which were selected as the best representative samples of each tissue type. HGCA2.0 presents subclades of coexpressed genes to a gene of interest, and performs various built-in gene term enrichment analyses on the coexpressed genes, including gene ontologies, biological pathways, protein families, and diseases, while also being unique in revealing enriched transcription factors driving coexpression. HGCA2.0 has been successful in identifying not only genes with ubiquitous expression patterns, but also tissue-specific genes. Benchmarking showed that HGCA2.0 belongs to the top performing coexpression webtools, as shown by STRING analysis. HGCA2.0 creates working hypotheses for the discovery of gene partners or common biological processes that can be experimentally validated. It offers a simple and intuitive website design and user interface, as well as an API endpoint.

A functional gene module identification algorithm in gene expression data based on genetic algorithm and gene ontology.

Since genes do not function individually, the gene module is considered an important tool for interpreting gene expression profiles. In order to consider both functional similarity and expression similarity in module identification, GMIGAGO, a functional Gene Module Identification algorithm based on Genetic Algorithm and Gene Ontology, was proposed in this work. GMIGAGO is an overlapping gene module identification algorithm, which mainly includes two stages: In the first stage (initial identification of gene modules), Improved Partitioning Around Medoids Based on Genetic Algorithm (PAM-GA) is used for the initial clustering on gene expression profiling, and traditional gene co-expression modules can be obtained. Only similarity of expression levels is considered at this stage. In the second stage (optimization of functional similarity within gene modules), Genetic Algorithm for Functional Similarity Optimization (FSO-GA) is used to optimize gene modules based on gene ontology, and functional similarity within gene modules can be improved. Without loss of generality, we compared GMIGAGO with state-of-the-art gene module identification methods on six gene expression datasets, and GMIGAGO identified the gene modules with the highest functional similarity (much higher than state-of-the-art algorithms). GMIGAGO was applied in BRCA, THCA, HNSC, COVID-19, Stem, and Radiation datasets, and it identified some interesting modules which performed important biological functions. The hub genes in these modules could be used as potential targets for diseases or radiation protection. In summary, GMIGAGO has excellent performance in mining molecular mechanisms, and it can also identify potential biomarkers for individual precision therapy.

Gene function and cell surface protein association analysis based on single-cell multiomics data.

Single-cell transcriptomics provides researchers with a powerful tool to resolve the transcriptome heterogeneity of individual cells. However, this method falls short in revealing cellular heterogeneity at the protein level. Previous single-cell multiomics studies have focused on data integration rather than exploiting the full potential of multiomics data. Here we introduce a new analysis framework, gene function and protein association (GFPA), that mines reliable associations between gene function and cell surface protein from single-cell multimodal data. Applying GFPA to human peripheral blood mononuclear cells (PBMCs), we observe an association of epithelial mesenchymal transition (EMT) with the CD99 protein in CD4 T cells, which is consistent with previous findings. Our results show that GFPA is reliable across multiple cell subtypes and PBMC samples. The GFPA python packages and detailed tutorials are freely available at https://github.com/studentiz/GFPA.